Fermionization and bosonization of expanding ID anyonic fluids 
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The momentum distribution of an expanding cloud of one-dimensional hard-core anyons is stud- 
ied by an exact numerical approach, and shown to become indistinguishable from that of a non- 
interacting spin-polarized Fermi gas for large enough times (dynamical fermionization). We also 
consider the expansion of one-dimensional anyons with strongly attractive short-range interactions 
suddenly released from a parabolic external potential, and find that momentum distribution ap- 
proaches that of its dual system, the ideal Bose gas (dynamical bosonization) . For both processes 
the characteristic time scales are identified, and the effect of the initial confinement is analyzed 
comparing the dynamics associated with both harmonic and hard-wall traps. 
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The recent realization of ultracold atoms in tight- 
waveguides has spurred extensive studies of one- 
dimensional (ID) systems exhibiting Bose-Fermi duality 
The inextricable link between interactions and ex- 
change statistics in ID allows to relate the Bose gas with 
(5-interactions with spin-aligned attractive Fermi gas with 
inverse coupling constant. The upshot is that, in spite of 
the different symmetry of the many-body wavefunction, 
dual systems share many of its properties, including any 
local correlation function. 

Paradigmatic examples are the ID gas of spin 
polarized-fermions and impenetrable bosons (the so- 
called Tonks-Girardeau (TG) gas 0) which have been 
prepared in the laboratory [3| . Similarly, the duality be- 
tween non-interacting bosons and spin-polarized fermions 
with strongly attractive odd- wave interactions (fermionic 
Tonks-Girardeau (FTG) gas [4]) has recently been ex- 
plored. A staggering manifestation of the BF duality is 
the dynamical fermionization in which the momentum 
distribution of the TG gas approaches during a ID free 
expansion that of non-interacting fermions [a, 0|- The 
reverse phenomenon has been predicted for FTG under 
the same dynamics, in which asymptotically the momen- 
tum distribution of the ideal Bose gas is exhibited [?[■ 
Far from being limited to the Bose-Fermi mappings, the 
duality in lower dimension extends to particles with frac- 
tional statistics, whose realization with ultracold atoms 
has recently been proposed [8] . Motivated by its applica- 
tion to transport in quantum Hall fluids, and the possibil- 
ity of performing universal quantum computation [3| , the 
properties of ID anyons are under current scrutiny. Par- 
ticularly, several studies have been devoted to ID anyons 
with short-r ang e delta interactions (anyonic Lieb-Linigcr 

and the limiting cases of im- 



model) i3,[iir[l2,[i3, 
penetrable anyons [15 



171 . |18|, though much less is 



known about its dynamics. 

In this letter, we shall consider the anyonic general- 
ization of TG and FTG gases, and show that they do 
undergo dynamical "dualization" during their free time 
evolution. 



Hard-core anyonic gas. Strongly (5-interacting 
anyons, at low enough densities may reach the so-called 
Tonks-Girardeau regime in which the particles become 
impenetrable. We assume that the gas of hard-core 
anyons (HCA) is initially confined in a harmonic trap of 
frequency ujq, whose eigenstates are denoted by (f>n{x). 
The description in such regime is simplified introducing 
a dual system, a spin polarized Fermi gas, whose 
ground-state wavefunction is a normalized Slater de- 
terminant '^Fixi,...,XN) = 7^det'j^~(^Q^j(/)„(x;). 
The wavefunction of the ATG gas is obtained from "ifp 
by imposing the correct symmetry under permutation 
of particle using the anyon-fermion (AF) mapping 15 1 
'^hca{xi, ■ ■ ■ ,xn) = AI(xi,. . . ,Xjy)-^p{xi,. ..,Xn), 



where Ae = Y\i<j^k<N one-parameter 
family of unitary operators which generalizes the 
well-known Bose-Fermi map 01 and e{x) = 1 (-1) if 
X > {< 0) and e(0) = 0. Varying the statistical 
parameter 9, the map smoothly extrapolates between 
spin-polarized fermions {9 = 0), and the well-known 
bosonic Tonks-Girardeau gas {9 = tt) [l^. Note that 
A^g — AI — Ag^ and the expectation values coin- 
cide {'i>HCA,'^HCA) = (*_f,*f), which imply that 
both dual systems share the same density profile [l3| 
pe{x,t) = N ^ \^{x,X2,. ■ ■ ,XN;t)\'^dx2- ■ -dxN = 
'^n=o \'t'n{x,t)\'^ ■ However, its momentum dis- 
tribution n{k) — {2Tr)^^ J J dxdye^'''^^~y^ p{x,y) 
is drastically different. For the spin polarized 
fermions, described by an 
wavefunction, the reduced 
sity matrix (RSPDM) simply 

N J '^f{x,X2, . . . ,XN)*'^Fiy,X2, ■ ■ 

'l2n=o 4>n{x)4>n{y) (normalized to N), whence it follows 
that npik) = ^^=0 l<^n(^)P i^^^ til^i^ denoting Fourier 
transform). For the HCA gas, making use explicitly 
of the anyon-Fermi map, the orthonormality of the 
single-particle eigenstates, and the Laplace expansion of 
the determinant in the dual wavefunction, it is found 



Slater determinant 
single-particle den- 
reads pF{x,y) — 
■ ,XN)dx2 - ■ -dxN = 
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FIG. 1; Momentum distribution of harmonically trapped 
hard-core anyons, smoothly extrapolating between the 
bosonic Tonks-Girardeau gas [6 = 0) and spin-polarized non- 
interacting fermions [6 — tt), for N — 5 particles (fco — 
y^mujo/h, in all figures). 



that the RSPDM can be efficiently computed as 



Phca 



l.n=0 



where A{x,y) — (P ^)-'"detP and the elements of the 
matrix P{x,y) are 

Pin = 6in - (1 - e-^'^y-'^My - x) r dz0r(z)0„(z) 



a result which holds under time-evolution. This gener- 
alizes for any anyonic parameter, the result recently ob- 
tained by Pezer and Buljan f20j for the Tonks-Girardeau 
gas (9 = tt). As shown in Fig. [1] an asymmetric mo- 
mentum distribution results for any statistical parame- 
ter other than = 0, tt fl3, [Hi- Moreover, for 9 = 0, 
Pin — Sin so the RSPDM becomes diagonal, and the 
familiar result for spin-polarized fermions is recovered. 
In what follows we shall be interested in the expansion 
dynamics after suddenly switching off the confining po- 
tential at i = 0. We find that, as the time of evolution 
goes by, the momentum distribution of the HCA gas ap- 
proaches that of its dual system, the spin-polarized ideal 
Fermi gas. 

Provided that Ae is time-independent (the quantum 
statistics is preserved under time-evolution), the expan- 
sion dynamics of the manybody wavefunction is found 
by using in the anyon-fermion mapping. The self-similar 
evolution of the single-particle states according to the 
time-dependent Schrodinger equation, j^^^nlSili 



dt 



2m dx^ 



'e(-t)- 



(pnix, t), with 8(t) being the Heav- 
iside step function, can be found exploiting the scaling 

law, ?i„(x,t) = 0„(x/5(t),O)e™^'''/2'''^-*-^"^(*)/'*/v/6(i) 
where the scaling factor b{t) is the solution of the differ- 
ential equation b+ijj'^{t)b = satisfying 6(0) = 1 and 
6(0) = 0,E„ = hujo{n+l/2), andT(<) = J*dt'/b^(t') 
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FIG. 2: Dynamical fermionization on the momentum distri- 
bution of hard-core anyons during a ID free-expansion for 
different statistical parameters (A^ = 5). Asymptotically, 
the momentum distribution of non-interacting spin-polarized 
fermions is obtained. The time of evolution is measured in 
units loq^ where wo is the frequency of the initial harmonic 
confinement. 



For uj{t) = a;o0(-i) it follows that b{t) = ^\ + L^y? , and 
therefore, for all t ^ loq^ the expansion becomes ballis- 
tic [23 | . During the expansion a dynamical fermion- 
ization occurs, see Fig. [21 by means of which the mo- 
mentum distribution tends to that of the non-interacting 
spin-polarized Fermi gas, np{k) = X]^=o^ l0n(fc)P where 
4>n{k) = (-l)"(2x^/^)i/4e-fc'-oi7„(v^fc^i7„ are the 
Hermite polynomials and xq = ^/h/mujo- Indeed, 
the stationary phase method (SPM) allows to find the 
asyrnptotic form of the momentum distribution for t — > 
oo [6], UAHAik) ~ \ujQ/b\nF{uJak/b) ^ np{k). Intu- 
itively, as the gas expands, the particles stop to in- 
teract, and asymptotically one find the distribution of 
quasi-momenta, which are integrals of motion. The 
fermionization time scale can be approximated astp (9) « 
sin(6'/2)/A^tJO: dependent on the statistical parameter 9 
and decreasing with the number of particles N. 

Expansion from a box. The scaling law is specific 
of the harmonic confinement and we may rightly won- 
der how the subsequent time-evolution differs from the 
expansion from other type of traps. In what follows 
we shall compare the previous results with the expan- 
sion from a box- like trap [23| . The main difference is 
that at variance with the parabolic potential the expan- 
sion from a hard- wall confinement is not self-similar, as 
the space-time dynamics already exhibits a rich tran- 
sient structure at the single-particle level. For a hard- 
wall confinement of width L, the well-known orthonor- 
mal eigenstates read </?j(a;) = {2/ Ly/^ sm{kjx)x[o.L]{x) 
{kj = jTr/L, j £ N, X[o,L]{x) is the characteris- 
tic function in the interval [0, L], and the correspond- 
ing energy eigenvalues are Ej = h^k'^/2m). The ex- 
act time-evolution of the single-particle eigenstates of 
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FIG. 3: Effect of the confining potential. The momentum 
distribution of the anyonic Tonks-Girardeau gas released from 
the box (solid lines) broadens with respect to the harmonic 
case (dashed lines), and undergoes fermionization in a shorter 
time scale (A'^ = 5). Note that uiq and L are chosen for the 
total energy of the gas to be the same in both traps (see text). 



a hard-wall trap in a ID free expansion (pj{x,t) can 
be found using the superposition principle Lpj{x,t) = 
J^^dx'Ko{x,t,x' ,t — 0)(pj{x',0) where the free propa- 
gator Xo(a;, t, a;', 0) = {m/2mhtY>'^exp{im{x~x'Y /2ht). 
After performing the integral explicitly, ipj{x,t) = 
Ea=±i"[e*"''''^M(a; - L,akj,T) - M{x,akj,T)] 
ht/m, the Moshinsky function is defined as 
1^ 

2 



iV2L 

where t 



M{x,k,t) = e*^u;[ - (fc - f ) ] /2 [2^, and the 

Faddeyeva function is related to the complementary er- 
ror function w{z) = e~^^Erfc(— iz) [2^. For all j > 1, the 
density profile \(pj{x, t)\'^ bifurcates in two main branches 
moving with ±fcj after a critical time tj = mL'^ /2jTrh re- 
flecting the underlying bimodal momentum distribution 
[2^ . For the free expansion of a TG gas the transition 
to the ballistic regime sharply takes place at In To 
compare the expansion from both harmonic and hard- 
wall traps we impose the particle number N and total 
energy of the AHC gas to be the same. For the many- 
body ground-state, the following relation must hold be- 
tween the frequency luq and the length of the box, L: 
LUo = hn'^iN + 1)(2N + l)/6mL'^N. Figure [3] shows how 
the momentum distribution is broadened with the initial 
hard- wall confinement with respect to the harmonic case, 
due to the dispersion relation En oc and the presence 
of higher momentum components. The reference case for 
6 = TT, the bosonic TG gas, presents a sharper peak at 
k = when prepared in the later trap. Note also how 
asymptotically the quasi-momentum distribution is al- 
ready flat (for = 5) for the box confinement, while in 
the harmonic trap case keeps an oscillating profile, both 
of which are eventually mapped into the coordinate den- 
sity profile in a, tp time scale. 

Anyonic attractive Tonks-Girardeau gas. The anyonic 
attractive Tonks-Girardeau (AATG) gas has recently 
been introduced in analogy with the Fermionic Tonks- 
Girardeau gas (FTG) The FTG results when a spin- 
polarized Fermi gas is driven through a p-wave Feschbach 
resonance inducing strongly attractive short range odd- 
wave interactions [27|. Indeed, in a given sector xi < X2, 
the resulting pseudopotential consists of a hard-core wall 
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FIG. 4; Momentum distribution of attractive hard-core 
anyons, smoothly extrapolating between the spin-polarized 
fermionic Tonks-Girardeau gas {6 = 0) and the non- 
interacting Bose gas {6 — n) for TV = 5 particles in a parabolic 
external potential. 



at X12 = xi — X2 = and a square well of depth V 
and width a in the limit a — > 0, — > 00, keeping V = 
[hiT / a)"^ / Am. The dual system of the FTG is simply the 
ideal Bose gas, whose ground state is simply described 
by a Hartree product ^ b{xi, . . . ,xn) = ]\n=i^oi^n), 
(j)o being the single particle ground state of the exter- 
nal potential. This motivates, the study of the ID 
anyons with such infinitely attractive pseudopotential - 
the AATG gas-, whose ground state wavefunction is given 
by the anyon-boson mapping ^/AATcixi, . ■ . ,xn) = 
ni<j<fc<Ar ^(^fc ~ ij)e"*2«(^fc-^3)5'3(a;i, . . .,xn), where 
the mapping holds under any unitary time evolution [l5| . 
A one-parameter family of RSPDM for the AATG results 
from direct computation in the line of [27] , 



PAATG{x,y) = N(t>o{x)(pQ{y)[PQ{x,y)] 



N-l 



(2) 



where 



P^{x,y) = / dze(y-z)e(a;-z)e-'i[^(^-^)-^(---)l|</.o(z)|^ 



1 - [1 + e-'''<y~''^] 



dz\Mz)\' 



(3) 



In particular, for the harmonic confinement, P^{x^y) 



— |Erf(y/xoV6) - Erf(a:/xoV6)|. The first 
term in Eq. N(l)Q{x)(f>o{y), corresponds to the 

RSPDM of the dual non-interacting Bose system. More- 
over, the momentum distribution ng(k) can be obtained 
by the usual double Fourier transform. Figure [4] shows 
ne(fc) for a gas of AATG in a harmonic trap as a function 
of the statistical parameter 9. The ideal gas distribu- 
tion N\(j)Q{k)\'^ is recovered from the general expression 
for 9 ~ TT, whereas the FTG gas is obtained at 9 = 0. 
However, for any other value of the statistical parameter 
9 the momentum distribution of the trapped system is 
asymmetric as a result of the anyonic permutation sym- 
metry. We shall next describe the one dimensional ex- 
pansion dynamics of a given AATG, following a sudden 
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FIG. 5: Dynamical bosonization on the momentum distribu- 
tion of attractive hard-core anyons during a ID free-expansion 
from a harmonic trap for different statistical parameters 
(A'^ = 5). Asymptotically the momentum distribution of an 
ideal Bose gas is recovered. 



switch-off of the hard- wall trap. Figure [5] shows how as 
the system expands in free space, the AATG gas un- 
dergoes bosonization reaching asymptotically the (sym- 
metric) momentum distribution of the ideal Bose gas, 
nsik) — N\(j)Q{k)\'^ . For t — s- oo the argument of each 
of the complementary error functions in P^{x,y) tends 
to vanish, and hence, Erf(O) = 0, and P^{x^y) — 1. It 
follows that PAATG{x,y) ~ N<pQ{x)(l)o{y), meaning that 
any one-particle observable, both local and non-local, 
becomes identical for both the AATG and ideal Bose 
gases. The SPM leads for the the asymptotic momen- 
tum distribution nAATcik) |wo/6|«s(ti'o^/&) ^ nsik). 
For NjUot 3> 1, the off-diagonal correlations vary expo- 
nentially as /O^atg (a;, y, i) ~ pB{x,y,t) exp[exp(-i6'e(y- 
a::)/2)iVcos(6'/2)|?/ — a;|/a;oiv^2;o)]. Hence, the bosoniza- 
tion time scale becomes tsiB) ~ N cos{9 /2)uJq^ , this is, 
tB{0) ~ N'^tp{TT — 6). Should we have considered the 
expansion from a hard-wall trap, the characteristic time 
scale would be tB(6') « NmL'^ cos{0/2)/2ttH. 

In conclusion, we have shown that ID hard-core anyons 
undergo dynamical fermionization during a free time- 
evolution, namely, the momentum distribution of the 
gas approaches that of spin-polarized non-interacting 
fermions. We have further shown that the complemen- 
tary dynamical process can take place: for anyons with 
strongly attractive short range interactions (the anyonic 
attractive Tonks-Girardeau gas), the freely time-evolving 
momentum distribution tends to that of an ideal Bose 
gas. In both cases the momentum distribution becomes 
symmetric, as the role of the anyonic permutation sym- 
metry diminishes. One may expect that the recently 
developed methods for the dynamics of a bosonic Lieb- 
Liniger gas 0|, could equally be generalized for frac- 
tional statistics allowing to study the expansion dynam- 
ics for finite interactions 1291 . 
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